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Genetic analysis of heterotic loci detected in a cross between indica and japonica 
rice {Oryza sativa L.) 
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The study on the genetic basis of heterosis has received significant attention in recent years. In this study, 
using a set of introgression lines (ILs) and corresponding testcross Fi populations, we investigated heterotic 
loci (HL) associated with six yield-related traits in both Oryza sativa L. subsp. indica and japonica. A total 
of 41 HL were detected on the basis of mid-parent heterosis values with single-point analysis. The F ; test- 
cross population showed superiority in most yield-related traits and was characterized by a high frequency 
of overdominant HL. Thirty-eight of the 41 HL were overdominant, and in the absence of epistasis, three HL 
were dominant, suggesting that heterotic effects at the single-locus level mainly appeared to be overdominant 
in rice. Twenty-four HL had a real positive effect, suggesting that they are viable candidates for the improve- 
ment of rice yield potential. Compared with the quantitative trait loci (QTLs) detected in the ILs, only six 
out of the 41 (14.6%) HL were detected in QTL analysis under the same statistical threshold, indicating that 
heterosis and trait performance may be conditioned by different sets of loci. 
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Introduction 

Use of heterosis has become a major strategy for increasing 
the productivity of plants and animals (Hua et al. 2003); 
however, the study of the genetic basis of heterosis has 
fallen far behind the exploitation of heterosis. Dominance 
(Bruce 1910, Jones 1917) and overdominance (Shull 1908) 
are two hypotheses proposed a century ago to explain the ge- 
netic basis of heterosis. Recent advances in genome research 
involving a number of molecular-marker techniques and the 
availability of high-density molecular linkage maps, togeth- 
er with developments in analytical methods (Lander and 
Botstein 1989, Zeng 1994), have facilitated the analysis of 
the genetic basis of quantitative traits. 

Recently, many quantitative trait loci (QTL) mapping 
studies have provided insight into the genetic basis of het- 
erosis (Frascaroli et al. 2007, Garcia et al. 2008, Hua et al. 
2003, Li et al. 2001, Li et al. 2008, Luo et al. 2001, Luo et 
al. 2009a, Mei et al. 2003, 2005, Melchinger et al. 2007, 
2008, Stuber et al. 1992, Xiao et al. 1995, Yu et al. 1997), 
and resulted in different explanations for this phenomenon, 
such as dominance, overdominance, epistasis and so on. We 
analyzed previous studies and found that differences in seg- 
regating populations and/or statistical methods caused the 
different conclusions. Furthermore, one known problem in 
establishing the genetic basis of heterosis has been the use of 
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whole-genome segregating populations, where interactions 
often mask the effects of individual loci (Semel et al. 2006). 

A solution to this problem is the use of introgression lines 
(ILs), in which small chromosomal segments are intro- 
gressed from the donor into the recurrent parent by consecu- 
tive backcrossing and selling (Eshed and Zamir 1994, 1995). 
Consequently, ILs provide more precise estimates of the ge- 
netic effects of introgression against a relatively uniform and 
elite lineage background (Tanksley and Nelson 1996) and 
are therefore well-suited for genetic analysis of heterosis. 
Analysis of overlapping chromosomal segments in ILs has 
proven to be a powerful strategy to more precisely map QTL 
and validate the QTL mapped in early generations or in 
genome-wide segregating populations (Paterson et al. 1990). 
In addition, several reports have indicated that ILs are a 
powerful tool for identifying new genes (Eshed and Zamir 
1994, 1995, He et al. 2006, Luo et al. 2009b), distinguishing 
pleiotropy from linkage (Yamamoto et al. 1998), eliminat- 
ing QTL linkage drag, and for map-based cloning (Alpert 
and Tanksley 1996). 

Recently, we analyzed heterotic loci (HL) between wild 
and cultivated rice associated with six yield-related traits in 
a set of 265 ILs, and our results supported the overdominant 
model involving a single functional Mendelian locus in the 
absence of epistasis (the results have been submitted to 
Genetics Research). In this study, we further investigated 
the genetic basis of rice heterosis using a set of 70 ILs from 
an intersubspecific cross. The lines were generated from a 
cross between IR24, a commercial indica cultivar, as the re- 
current parent, and Asominori, a typical japonica cultivar, as 
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the donor parent. On the basis of the set of 70 ILs and 175 
testcross F|S (derived from the cross between the ILs), QTL 
and HL associated with yield and yield-related traits be- 
tween the ILs and testcross FiS were analyzed. The genetic 
effects and main features of the HL are discussed. 

Materials and Methods 



markers were chosen to definitively identify the F] hybrid 
plants. Five plants from the middle of each row were har- 
vested at maturity in both the ILs and F ; individually, and 
the following traits were scored: spikelet number per panicle 
(SP), filled grain number per panicle (GP), percent seed set 
(SSP), 1000-grain weight (GW), panicle number per plant 
(PP) and grain yield per plant (YP). 



Experimental population development 

In this study, an IL population composed of 70 lines car- 
rying variant introgressed segments of Asominori (japonica) 
was used. The IL population included the indica cultivar 
IR24 background and covered the entire japonica genome 
with overlapping introgressed segments of each line. The de- 
tailed characteristics of ILs were presented in Aida et al. 
(1997) and Kubo et al. (2002). Only the heterozygous intro- 
gressed segments were all eliminated through associating 
ILs selfing and corresponding markers analysis. Additional- 
ly, 175 F i testcross individuals, derived from crosses among 
the 70 ILs as described below, were evaluated. 

Crosses were made between lines chosen by random per- 
mutations of the 70 ILs. In each round of permutation, the 70 
ILs were randomly divided into two groups, and the lines in 
the two groups were paired up at random to provide the par- 
ents for 35 crosses. Each of the 70 ILs was used only once in 
each round of pairing and crossing. This procedure was re- 
peated five times, resulting in a population consisting of 175 
crosses. The design resembled that previously reported by 
Hua et al. (2003). Crosses to generate the F ; generation were 
completed in the summer of 2004 in Beijing and the winter of 
2004 in Sanya (18°20'N, 109°50'E, Hainan Province, China). 

Field trials and phenotypic evaluations 

Two separate experiments were conducted at the Experi- 
ment Stations of China Agricultural University (CAU, 
Beijing, 39°54'N, 116°28'E) and the Anhui Academy of 
Agricultural Sciences (AAAS, Hefei, 3P51'N, 1 17°18'E) in 
2005. In the CAU environments, the ILs, parents, and test- 
cross FiS were sown in the seedling nursery on May 7, 2005. 
The 25-d-old seedlings were transplanted into 6-row plots 
(72 plants), each consisting of two rows of Fj hybrids and 
the corresponding two IL parents in fixed order: IL ; , Fi and 
ILj (j = I + n, namely j > I, n = 1, 2, , 69). Planting of ex- 
perimental units followed a complete randomized block de- 
sign with two replications. All materials were planted with 

13.3 cm spacing between plants within each row and 

26.4 cm spacing between rows. Field management followed 
normal rice production conditions. In the AAAS environ- 
ment, all of the materials were sown on June 5, 2005 and 
transplanted into the field after 25 d. The field arrangement 
was the same as the CAU environment. 

Morphological features characteristic of hybrid rice were 
used to identify the Fj hybrids, including plant height, head- 
ing date, and tillers, among others. Parents were used as the 
baseline for comparison, and observations were made 
throughout the growing season. Several appropriate RFLP 



DNA extraction and molecular marker analysis 

DNA was extracted from fresh leaves according to the 
CTAB method (Murray and Thompson 1980) with minor 
modifications. The RFLP markers analyzed in this study 
were from a previous publication (Saito et al. 1991). The 
molecular marker order was based on the rice linkage map 
described by Saito et al. (1991) and Gramene (http://www. 
gramene.org/db/markers/marker_view) (Fig. 1). A total of 
87 polymorphic RFLP markers were used to genotype the 70 
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Fig. 1. Genomic locations of putative HL associated with yield- 
related traits detected using single-point analysis in the testcross F^. 
Underlined markers correspond to HL at right. 
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ILs and the recurrent parent IR24. The F ; testcross geno- 
types were deduced on the basis of the genotype of its corre- 
sponding parental IL. 

Analysis method of mapping HL 

The direct trait measurement values from the six yield- 
related traits obtained from the ILs were used to identify the 
associated QTLs. The mid-parental heterosis values, H M? = 
Fi - (IL; + ILj)/2, (where F, is the testcross Fj mean value for 
the measured trait, and IL ; and ILj are the corresponding IL 
parent mean values for the same measured trait) were used to 
identify the loci affecting heterosis (HL is a locus that dem- 
onstrates a significant difference between the heterozygote 
and the mean of two corresponding homozygotes; i.e., HL is 
the QTL for heterosis) in the six yield-related traits. The test- 
cross F] trait measurements were used to identify the loci af- 
fecting testcross Fj performance. 

On the basis of the IL structure, QTLs or HL can be 
mapped on introgressed chromosome segments. One repre- 
sentative marker of each specific introgressed segment was 
defined as a QTL or HL. The association between the pheno- 
typic and RFLP data was investigated by single-point analy- 
sis using the software package Map Manager QTXM7 
(Manly et al. 2001). 

The genetic effects of HL from IL populations (advanced 
backcross population) were analyzed by single-locus associ- 
ation using QTXM7. The statistical a priori threshold for 
main effect loci was P < 0.01. If a group of similar markers 
was associated with the same trait of similar effect(s) in both 
magnitude and direction, the marker (or locus) exhibiting the 
lowest P-value was chosen as the locus associated with the 
trait(s). 

The genetic effects in the testcross Fjs were defined as 
follows: d = H MP = [Fi - (ILj + ILj)/2]; the trait mean values 



in the testcross F t s were Fj = (a + d), where a represents the 
additive effects from the performance values of testcross F|. 
Subsequently, HL effects were inferred by comparison of 
the genetic effects on Fi performance and mid-parent hetero- 
sis. HL with d/a < 1 were considered complete or partially 
dominant loci, and expected to generate an estimate of Fi 
performance (a + d) equal to or higher than twice the //mp 
(d). HL with d/a > 1 that is, 2d (2 x H M? ) >a + d(F{), or only 
detectable for H MP were determined overdominant loci (Li et 
al. 2001, Mei et al. 2005, Melchinger et al. 1998). 

Results 

IL and Fj testcross performance and Hup for yield-related 
traits 

The phenotypic data from the ILs and Fi testcrosses 
showed a continuous distribution (Supplemental Fig. 1), 
which suggested that the yield-component traits were quan- 
titative traits under multi-gene control. The performance sta- 
tistics for ILs and testcross F[S, and H M? of Fi testcrosses 
between two environments are provided in Table 1 . Com- 
pared with ILs, testcross FjS had a significantly (P<0.05) 
greater mean value for each trait in the environment of the 
China Agricultural University (CAU) Experiment Station 
(Table 1), and clearly demonstrated hybrid vigor. In the envi- 
ronment of the AAAS Experiment Station, the F| mean value 
for each trait was also greater than that of ILs, while the dif- 
ferences were not significant (P < 0.05) for the GW and SSP 
traits. The highest heterosis was detected in grain yield per 
plant, while 1000-grain weight showed the lowest heterosis 
in both environments, consistent with the results from other 
studies (Hua etal. 2003, Li etal. 2001, Luo et al. 2001). The 
heterosis levels for all yield-related traits varied in the 175 
testcross F,s, from highly negative to highly positive. 



Table 1. Summary of performance statistics for ILs, testcross FiS and H MP of F| testcrosses between two environments 



Locations 


Trait" - 


Performance of ILs 


Performance of testcross F t s 






of testcross F]S 






MeantSD 


Range 


Mean±SD 


Range 


Mean±SD 


SE 


% c 


Range 


% d 


Beijing 


SP 


92.9 ±17.8 


52.7-128.7 


104.3 ± 16.0* 


66.2-168.8 


10.9± 12.0 


0.91 


3.0 


-7.7-100.2 


-1.9-36.5 




GP 


76.4 ±14.4 


44.7-108.8 


92.1 ±14.8* 


56.7-153.2 


14.3 ±11.5 


0.87 


4.7 


-6.4-85.4 


-2.0-35.7 




SSP (%) 


82.6 ±7.4 


54.7-94.4 


88.3 ±4.8* 


58.2-96.0 


4.4 ±4.9 


0.37 


1.5 


-12.8-20.3 


-4.5-7.6 




GW(g) 


22.9 ±2.9 


16.3-33.2 


23.9 ±1.7* 


19.0-28.3 


1.0± 1.4 


0.10 


1.1 


-3.7-5.8 


-3.6-7.4 




PP 


10.0± 1.5 


6.3-13.1 


11.3 ±2.0* 


3.2-21.0 


1.3±2.0 


0.15 


3.4 


-6.2-10.7 


-16.4-25.7 




YP (g) 


18.1 ±3.3 


10.7-24.9 


25.1 ±5.7* 


8.9-48.0 


6.8±5.9 


0.45 


9.9 


-10.3-30.0 


-13.4-43.8 


Hefei 


SP 


121.0 ± 15.7 


86.5-160.4 


130.2± 17.3* 


82.9-188.1 


9.5 ±13.4 


1.00 


2.0 


-15.8-63.4 


-3.7-12.7 




GP 


96.5 ± 16.4 


45.0-131.9 


104.3 ± 14.7* 


70.4-146.4 


8.1 ±13.3 


1.00 


2.3 


-35.6-41.3 


-8.4-12.8 




SSP (%) 


79.8 ±9.5 


50.3-95.9 


80.4 ±7.2 


58.4-93.3 


0.9 ±8.5 


0.64 


0.3 


-18.0-23.7 


-7.2-9.4 




GW(g) 


22.8 ±2.1 


18.3-28.5 


22.5 ± 1.4 


19.0-26.9 


0.2 ± 1.0 


0.54 


0.1 


-2.2-2.4 


-2.1-2.7 




PP 


1 1.0+ 1.3 


8.7-14.5 


13.5±2.5* 


8.4-28.0 


2.4 ±2.6 


0.07 


5.6 


-3.1-17.0 


-6.7-38.6 




YP (g) 


25.3 ±4.5 


12.5-37.0 


33.3 ±7.2* 


18.1-64.6 


7.8± 7.1 


0.19 


7.9 


-6.8-38.2 


-6.2-41.9 



" Trait abbreviations: spikelet number per panicle (SP), filled grain number per panicle (GP), percent seed set (SSP), 1000-grain weight (GW), 

panicle number per plant (PP), grain yield per plant (YP). 
4 Mid-parental heterosis, H Me - Fj - MP, where MP was the mid-parental trait value (ILj + ILj)/2, Fi was the performance value of testcross Fi. 
c Rate of mid-parental heterosis = (F, - MP)/MP*100%. 
d Range of mid-parental heterosis rate. 

* P < 0.05 significantly different between the mean performance of testcross F[S and ILs. 
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Table 2. Phenotypic correlation (R) for six yield- related traits be- 
tween IL and testcross Fi performance values and// MP values between 
two environments 



Location Trait" 



Performance 
values of IL 6 
and testcross F[ 



H M? values and 
performance 
values of 
testcross F, 



Performance 
values of IL 4 
and H MP values 
of testcross F, 



Beijing 


SP 


0.663 


0.675 


-0.104 




GP 


0.636 


0.695 


-0.113 




SSP 


0.470 


0.502 


-0.527 




GW 


0.670 


0.376 


-0.436 




PP 


0.238 


0.866 


-0.280 




YP 


0.135 


0.902 


-0.306 




Mean 


0.469 


0.669 


-0.294 


Hefei 


SP 


0.637 


0.730 


-0.061 




GP 


0.508 


0.663 


-0.307 




SSP 


0.353 


0.693 


-0.430 




GW 


0.721 


0.374 


-0.373 




PP 


0.127 


0.923 


-0.264 




YP 


0.259 


0.909 


-0.169 




Mean 


0.434 


0.715 


-0.267 



" See Table 1 for abbreviations. 

* Performance values of ILs are indicated average performance values 
of IL; and ILj. 



Relationships among the mean trait values of ILs, Hmp, and 
Fj performance 

Table 2 shows the correlation coefficients among test- 
cross Fi mean values, mid-parent heterosis values, and pa- 
rental IL mean values for yield-related traits. Highly positive 
correlations between H M? and Fj performance were found 
for most traits in both environments; the average R (correla- 
tion coefficient) was 0.669 in CAU environments and 0.715 
in AAAS environments. A general negative correlation 
trend was evident between IL and H MP trait values; the aver- 
age R was -0.294 in CAU environments and -0.267 in 
AAAS environments. A positive correlation between IL trait 
values and the F[S was observed, while the correlation co- 
efficients appeared to show a larger difference among traits, 
with an average R of 0.469 (range of 0.135 for YP to 0.670 
for GW) in CAU environments and 0.434 (range of 0. 127 for 
PP to 0.721 for GW) in AAAS environments. 

QTLs for 6 yield-related traits in ILs performance, F/ test- 
cross performance, and heterosis (Hmp) 

The IL phenotypic data for the six yield-related traits 
were used to identify the associated QTLs. Twenty-five 
QTLs were detected (Table 3), but only two of the 25 QTLs 
were common between the two environments. Trait pheno- 
typic values from F| testcrosses and estimated H MP values 
were used to infer the QTLs contributing to F ; testcross per- 
formance and heterosis. Forty-eight QTLs influencing Fi 
testcross performance and 41 HL associated with// M p values 
were detected for the yield-related traits (Table 3 and Fig. 1). 
Nine of 48 F, performance QTLs and two of the 41 HL were 
common to the two environments. 



QTLs for spikelets per panicle: Six QTLs were identified 
in ILs in the two environments, with two QTLs were identi- 
fied in both environments, and three QTLs contributed to an 
increasing effect. Five HL and 11 loci affecting F ; perfor- 
mance were detected in testcross FiS, and four of these loci 
(near markers CI 12, C246, CI 172 and X249) were simulta- 
neously detected in both cases. Of these loci associated with 
F| performance, only three (near markers X249, R1877 and 
C1350) were identified in both environments and seven loci 
appeared to have positive effects on SP. Interestingly, all 
five HL exhibited a positive effect, and only one locus (near 
marker CI 12 on chromosome 1) was also detected by QTL 
analysis of ILs. Furthermore, a comparison of the genetic ef- 
fects of loci detected in both Hmp and Fj testcross perfor- 
mance indicated a d/a > 1 in qSpIlb and qSp3b\ suggesting 
overdominant loci and dominant effects in qSpl and qSp5. 
The qSp5' was only associated with mid-parent heterosis and 
showed overdominant expression. 

QTLs for filled grains per panicle: Four QTLs were de- 
tected in ILs in the two environments and two enhanced the 
phenotype. Seven HL and 10 loci affecting F ; performance 
were detected in testcross F|S and four were detected in both 
cases. Of these loci associated with mid-parent heterosis 
and/or Fj performance, nine loci had significant positive ef- 
fects and only four loci appeared to have negative effects. 
Furthermore, one locus influencing both mid-parent hetero- 
sis and Fi performance, near marker CI 12 on chromosome 
1 , was also detected in QTL analysis of ILs, and contributed 
an increasing effect in ILs and positive heterosis in testcross 
FiS on filled grains per panicle. Notably, all the HL appeared 
overdominant. 

QTLs for percent seed set: Three QTLs were detected in 
ILs in both environments, and only one derived from 
japonica showed a positive effect on seed set percent. Eight 
HL and five loci affecting F| performance were detected in 
testcross FiS, and only one was detected in both cases and 
appeared to be overdominant. The remaining seven HL were 
only mapped in mid-parent heterosis and exhibited overdom- 
inance. Only one HL was common to the two environments, 
and only one overlap between HL and IL QTLs was detected. 

QTLs for 1000-grain weight: Five QTLs were detected in 
ILs in both environments, and only one exhibited a positive 
effect. Ten HL and 12 QTL affecting Fj performance were 
detected in testcross Fjs and two were detected in both cases. 
Furthermore, nine out of 10 HL appeared to be overdomi- 
nant and the remaining locus (qGw3b) appeared dominant. 
Interestingly, the HL, qGwlb, was also detected in QTL 
analysis of ILs and showed a decreasing effect when homo- 
zygous, while showing a significant positive heterosis effect 
when heterozygous. 

QTLs for panicles per plant: In total, four QTLs were de- 
tected in ILs in the two environments, and two contributed a 
positive effect. Seven HL and seven QTLs affecting Fi per- 
formance were mapped in the testcross population, and three 
were detected in both cases. Furthermore, all 7 HL appeared 
to be overdominant, and six showed a significant positive 
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Table3. QTLs affecting six yield-related traits identified in ILs, testcross FiS and H M? 



Traits" Locations 


QTLs 


Markers - 




IL 






F, 










ji ii 
- dlcfi 




PV C 


pc 


a c 


PV C 


pc 


a + cf 


PV C 


pc 






i>F Beijing 


qSpl 


CI 12 


8 


0.0095 


10.1 1 


8 


0.0001 


15.96 


4 


0.0066 


8.02 


1.01 




qSp2 


Till m 








4 


0.0065 


8.07 












qbp3a 


CD 15 


Id 


O AO 1 £ 

0.U01D 


1 O TO 

— ly.zy 


















qSp3b 


X249 








5 


0.0021 


8.99 












qSp4 


R2373 








4 


a Mil 
0.00 11 


1 1 .44 












qSp5 


C246 








10 


0.0000 


34.98 


5 


0.0021 


17.45 


1.00 




qSp9 


CH63 








4 


0.0051 


-8.85 












qSplO 


R1877 








8 


0.0001 


—12.53 












qSplla 


C1350 








4 


0.0070 


-11.49 












qSpl lb 


CI 172 








6 


A AAAT 

U.OOU/ 


24.22 


3 


0.0091 


13.55 


1.12 


Hefei 


qSpla ' 


X92 


12 


0.0043 


1 1.13 


















qSplb ' 


CI 12 


9 


0.0096 


9.80 


















qSp3a ' 


C515 


11 


0.0057 


-15.02 


















qSp3b ' 


X249 








4 


0.0096 


8.89 


10 


0.0001 


1 1.22 


2.52 




qSp5 


X3S / 














1 A 
1U 


A AAA 1 
0.UUU1 


1 A C7 
1U.D / 






qSp6' 


X342 


1 1 


0.0075 


9.45 


7 


A AAAA 

u.oouy 


12.93 












qopa 


X278 








o 
5 


A AAAO 
U.UUUS 


— 2o.2U 












qSplO' 


R1877 


13 


0.0031 


-13.16 


10 


0.0001 


-16.40 












qopl 1 


C135U 


•7 


o oooo 

u.uuyu 


o <o 

— y.DO 


1 £ 

ID 


A AAAA 
U.OOUO 


—25.36 










CjF Beijing 


qGpl 


CI 12 


8 


A AAAO 

0.U098 


8.04 


8 


0.0001 


13.96 


7 


A AAA"} 

0.UUU3 


10.12 


1.45 




qGp3a 


C515 


15 


0.0014 


-15.81 


















qGp3b 


X249 








8 


A AAAO 

U.0002 


10.12 


6 


0.001 1 


6.86 


1.36 




qGp4a 


K23 /3 








4 


A AA*7A 

U.00/9 


10.47 












qGp4o 


POO 1 

c»yi 


1 A 

14 


A AAT) 

U.UU22 


o o/: 

y.»6 


















qGp5 


C246 








1 1 


A AAAA 

U.OOUO 


32.60 


6 


0.0010 


18.82 


1.15 




qGplv 


X3 / 








A 

y 


A AAA 1 

U.0UU1 


— 1 1.1D 












qGpll 


C1172 








5 


U.0015 


20.92 


3 


0.0094 


12.88 


1.23 




qGpl2 


X402 














4 


0.0056 


-6.66 




Hefei 


qGpl ' 


X92 








5 


0.0056 


1 1 .23 












qGp3 ' 


X249 








4 


0.0099 


6.47 












qGp5' 


X387 


8 


0.0096 


-8.43 








5 


0.0091 


7.27 






qGp6a ' 


X209 














5 


A flflTl 

O.UU/2 


-7.85 






qGpoo 


X342 








D 


A AAC 1 

U.OUD 1 


O A A 

s.4y 












qGpti 


CI 121 








o 

0 


A AAAC 
U.OUUD 


15. ID 












qGplO' 


X37 








7 


A AAAA 

u.uouy 


-9.99 












qGpl 1 


pi 
C3 








6 


a finn 
U.0U22 


1 A A 1 

— 14.43 










CCD [5 

sar Beijing 


qSspl 


vm 

xyz 


1 o 

lo 


A AAAC 
0.U00D 


—6.32 


2 


A AAA"7 

u.uuy / 


—2.36 












qSsp2a 


G227 








12 


A AAAA 

u.uouu 


-5.05 












qSsp2b 


X45 


10 


0.0096 


—4.28 








5 


0.0021 


3.1 1 






qiisp3a 


K3156 








1 o 

lo 


A AAAA 

u.uuuu 


nil 
— 12. 1 1 












qSsp3b 


X249 














5 


0.0025 


2.73 






qSsp4 


POO 1 

c»yi 














4 


a aa/;o 

o.uu&y 


O A A 

—2.44 






qSsp6 


X27 














2 


A AAAA 

u.uuyu 


-1.31 






qSsp? 


VI "2 O 

X338 














A 

4 


A AAC/i 

U.UU56 


3. /5 






qbsplu 


R1877 














4 


0.0096 


-2.53 






qSsp 1 2 


X402 














4 


0.0051 


-2.88 




Hefei 


qSspS ' 


C246 








7 


0.0018 


-15.70 


7 


0.0016 


-16.40 


2.09 




qSspSa ' 


X209 


12 


0.0044 


6.02 


















qSsp6b ' 


X27 














5 


0.0045 


-3.24 






qSspll ' 


C1172 








5 


0.0051 


-9.03 










GW Beijing 


qGwla 


C955 














8 


0.0001 


1.81 






qGwlb 


X92 


21 


0.0001 


-2.68 








18 


0.0000 


1.93 






qGw2 


C560 








4 


0.0068 


1.11 












qGw3a 


R19 








10 


0.0000 


-1.47 












qGw3b 


X249 








10 


0.0000 


-1.30 


3 


0.0094 


-0.59 


0.91 
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Table 3. (continued) 



Traits'" Locations 


QTLs 


Markers - 




IL 




F, 










- a/a" 




PV C 


pc 


a c PV C 


pc 


a + d° 


PV C 


pc 


a* 






q&w4 
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A 
4 


A AA^A 


1 AO 

— i.uy 


o 
o 


A AAA1 


1 T> 
— 1 .ZZ 


2.24 




qhwj 


A3 5 / 






4 


0.00o4 


A OA 

— o.yo 












qGw6a 


AZliy 






3 


0.0090 


A OA 
0.80 












qGwob 


V5 A O 

A34z 












1 c 
ID 


A AAAA 
0.0000 


1 A A 

1.44 






qGw8a 


CI 121 


12 


0.0058 


—2.29 5 


0.0035 


— 1.34 












q(jW6D 


X56 






9 


A AAAA 

o.oooo 


-1.71 












qGwlO 


R1877 






4 


0.0072 


0.89 












qGwlz 


cyot 






4 


A AAO 1 

o.ooyi 


A Cll 

— o.y3 










Hefei 


qGwla ' 


R210 






5 


0.0076 


0.89 












qGwlb ' 


Ayz 












o 


A AAO 1 

0.0021 


a n i 
0.91 






qGw3a ' 


C515 












6 


0.0034 


-1.01 






qGw3b ' 


R19 






6 


0.0027 


-0.94 












qGw5 ' 


A3 5 / 






0 


A AAOA 

0.0029 


A OO 

— o.sy 












qGw6' 


X209 


10 


A AAO O 

0.0088 


1.25 8 


A AAAO 

0.0008 


1.05 












qGw8a ' 


CI 121 












1 


A AA 1 C 

0.0015 


1.21 






qGw8b ' 


X56 


1 1 


0.0067 


-1.64 3 


A AAA"> 

0.0092 


-1.05 












qGwlla ' 


AJzU 












1 1 
1 1 


A AAAA 
0.0000 


1.24 






qGwllb ' 


CI 172 


9 


0.0097 


-1.86 9 


0.0002 


-2.35 












qGwlza 


Alov-z 












5 


0.00o4 


A OA 

0.99 






quwlZb 


rod i 












/: 
O 


A AA1 1 
U.UUJ 1 


A HZ. 

—0. /D 




PP Beijing 


qPpl 


CI /» 






J 


A AA1 C 

0.0035 


— 1.1 / 












qPp2a 


X45 


17 


0.0007 


— 1.13 
















qPp2b 


C560 






5 


A A AO A 

0.0030 


— 1.46 












qPp3 


X249 






3 


A AAAO 

0.0098 


-0.92 


5 


0.0017 


-1.16 


2.52 




qPp7 


A33o 






3 


A AAIO 

0.003v 


1 £ o 

i.5y 












qrplua 


A3 / 






5 


A AATA 
0.0020 


1. 10 


5 


A AAO 1 

O.OOzl 


1 n 
1.1 / 


") AO 

2.02 




qrplllb 


csoy 












A 

4 


A AACA 
0.0050 


A CI 






qPpluc 


C405 






5 


0.0034 


0.88 












qPplla 


X320 


1 1 


0.0079 


1.14 
















„ r> 17 1- 

qrpl lb 


V 1 OO 1 

Alov-I 












4 


A AAO/t 

o.ooy4 


1 OA 

1.29 




Hefei 


qPpS' 


R1553 






8 


A AAA/1 

0.0004 


2.83 


7 


0.0012 


2.67 


1.89 




qPpT 


X338 


1 1 


0.0055 


1.06 
















qPplO' 


R1629 












4 


0.0099 


1.42 






qPplV 


X52 












10 


0.0001 


6.98 






qPpl2' 


X402 


15 


0.0014 


-0.96 














YP Beijing 


qYpl 


C112 






3 


0.0092 


3.04 


6 


0.0006 


4.91 


3.23 




qYp4 


C1016 












5 


0.0037 


-6.09 






qYpl 


RH45 


12 


0.0050 


3.23 
















qYpl2 


R3375 






7 


0.0003 


3.15 










Hefei 


qYpl' 


R210 


11 


0.0070 


3.02 
















qYp6' 


X209 


24 


0.0000 


4.04 






4 


0.0091 


-3.77 






qYpll ' 


X52 






6 


0.0026 


15.25 


7 


0.0009 


16.41 


2.15 



" See Table 1 for abbreviations. 

* H M? is the mid-parental heterosis of testcross F, calculated from H M? = F[ - MP, where MP = (ILj + ILj)/2. 

0 PV, phenotypic variance explained by the locus; P, probability that the marker genotype had no effect on the trait; a, additive effects from the per- 
formance values of ILs; a + d, additive and dominance effects from the performance values of testcross Fi; d, dominance effect from H M? values. 

d d/a, ratio of dominant and additive effects. HL with d/a < 1 were referred to as complete or partial dominant loci. By contrast, HL with d/a > 1 
were referred to as overdominant loci (Melchinger et al. 1998, Li et al. 2001, Mei et al. 2005). 



heterotic effect on panicles per plant. Interestingly, no HL or 
F] performance QTLs were detected in the QTL analysis 
above. 

QTLs for grain yield per plant: Three QTLs were detect- 
ed in ILs in the two environments, all of which contributed 



increasing effects. Four HL and three F ; performance QTLs 
were detected in testcross Fis and two were detected in both 
cases. All of the HL appeared to be overdominant. No HL or 
Fi performance QTLs were detected in the QTL analysis 
above. 
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Table 3 indicates the genetic overlap of H MP and QTLs 
detected in the F ; testcrosses. Analysis detected 48 QTLs in- 
fluencing Fi performance; 16 were associated with H M? and 
32 QTLs with additive and dominant effects were detected 
for Fi performance. At the single locus level, 38 of the 41 
HL (92.7%) were overdominant and three appeared domi- 
nant. In IL QTL analysis, six of the 41 HL (14.6%) were re- 
solved at the same statistical threshold and showed less 
genetic overlap with the yield-related trait QTLs. 

Discussion 

The complex nature of heterosis makes it difficult to divide 
into individual components, particularly in F 2 , recombinant 
inbred and backcrossed populations, largely because of epi- 
static interactions among the many segregating loci through- 
out the genome (Li et al. 2001, Luo et al. 2001, Semel et al. 
2006); therefore, it is difficult to define specific heterotic 
phenotypes and the individual genomic loci that control 
them. However, the IL population allowed us to partition 
heterosis into defined genomic regions, largely eliminating 
genome-wide epistasis. To understand the action types of 
heterotic loci, QTLs and HL associated with yield and yield- 
related traits were investigated using a set of 70 japonica 
introgression lines against an indica background and 175 
testcross FiS derived from the crosses between the ILs. A 
total of 41 HL for six yield-related traits were identified. The 
heterotic effects were determined as the combined effects of 
both additive and dominant gene actions, estimated from the 
performance values of testcross Fjs and the dominance ef- 
fects estimated from their mid-parent heterosis (Hmp) values. 
On the basis of this strategy, we characterized the gene ac- 
tion type at each HL, and the 41 HL revealed two different 
genetic effects, dominant or overdominant. These HL data in- 
dicated that overdominance was the major underlying factor 
of heterosis in the absence of epistasis. Notably, Semel et al. 
(2006) carried out quantitative genetic and phenotypic anal- 
yses on an IL population of tomato (Solatium lycopersicum) 
carrying a single chromosomal segment from the distantly 
related wild species Solarium pennellii. In the absence of 
epistasis, at a single locus level, overdominant loci had 
greater effects on tomato yield and fitness. 

To analyze a set of 265 ILs and their testcross Fi of 
O. rufipogon Griff, against the background of the indica 
high-yielding cultivar Guichao 2 (O. sativa L.), Luo et al. 
(2011) identified 42 HL associated with six yield-related 
traits in wild and cultivated rice. Furthermore, 38 were over- 
dominant and, in the absence of epistasis, four HL were 
dominant (the results have been submitted to Genetics 
Research). To confirm the validity of overdominant HL, 
Luo (2006), Luo et al. (2011) evaluated several key HL, 
such as hyp2 (near RM236 on chromosome 2) and hspll 
(near RM224 on chromosome 11). The two loci coincided 
with the locations of qGY2-l and qGYll-2, respectively, 
two yield-improved QTLs mapped by Li et al. (2002). Luo 
(2006) and He et al. (2006) narrowed qGY2-l (correspond- 



ing to hyp2) down to a 102.9-kb region by constructing NILs 
that differed in only a single QTL. Notably, Luo (2006) con- 
firmed that qGY2-l was an overdominant heterotic locus. 
At another locus (qGYll-2), Luo et al. (2009b, 201 1) identi- 
fied a QTL and/or HL (hspll) associated with the number of 
spikelets per panicle, and narrowed hspll down to a 702-kb 
interval between markers RM224 and RM3577 using chro- 
mosome fragment substitution analysis. Genetic analysis 
indicated that hspll appeared to be weak overdominant ef- 
fects. In this study, the reproducibility of key lines and HL 
need to be further evaluated. Investigating F 2 progeny de- 
rived from the Fi of target ILs to check whether there was 
epistasis between introgressed regions and confirming the 
genetic effect of heterotic loci will be our next research. 
Therefore, our present study favored, in the absence of 
epistasis, at a single locus level, overdominance as major ge- 
netic basis of HL in rice. To confirm this conclusion, further 
work, such as fine mapping and cloning, needs to be per- 
formed. 

Many QTL mapping studies have shown that genetic ef- 
fects in Fi hybrid rice appear to fit one of three different 
models: additive, dominance and overdominance at a single 
locus level (Hua et al. 2003, Li et al. 2001, Li et al. 2008, 
Luo et al. 2001, Luo et al. 2009a, Mei et al. 2003, 2005). 
Our results are consistent with previous studies. We detected 
48 QTLs influencing F] performance, and 16 were associat- 
ed with H UP and revealed dominant or overdominant effects, 
the remaining 32 appearing to have additive and dominant 
effects. To analyze the H MP data, we identified 41 HL. Of 
them, 25 loci only influenced H MP and appeared to be over- 
dominant. The cause of heterosis of intersubspecific hybrids 
or intrasubspecific hybrids might show certain differences. 
In intersubspecific hybrids, many researches have shown 
that the overdominant effects were more advantageous than 
dominant effects (Li et al. 2001, Luo et al. 2001, Mei et al. 
2003, 2005). Li et al. (2008) simultaneously analyzed the 
heterosis of intersubspecific (IJ) hybrids and intrasubspecif- 
ic (II) hybrids using backcross populations of RIL, and 
found that the ratio of overdominant loci in II hybrids was 
higher than in I J hybrids. 

Li et al. (2001) and Luo et al. (2001) analyzed the corre- 
lation between RILs and backcrossed populations for yield- 
related traits, and considered that backcross Fi performance 
was mainly determined by dominant gene action. A sub- 
sequent study (Mei et al. 2005) supported evidence that 
backcross Fj performance was largely determined by non- 
additive gene action. Our study employed a different exper- 
imental design; random crossing between ILs to establish a 
testcross population. We investigated the correlation coeffi- 
cients among testcross F] mean values, mid-parent heterosis 
values and parental IL mean values, and found a highly pos- 
itive correlation between testcross F ; performance and H M p 
values and a lower positive correlation between IL and Fj 
testcross performance values for PP and YP in both environ- 
ments (Table 2). This indicated that dominant gene action 
rather than additive gene action was a substantial contributor 
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to Fi testcross performance for traits PP and YP, and these 
results were consistent with the high H MP values in the test- 
cross Fis. In contrast, a lower positive correlation between 
testcross F] performance and H MP values and a highly posi- 
tive correlation between IL performance values and F] test- 
cross performance values were found with the GW trait, 
suggesting that additive gene action was a substantial con- 
tributor to Fi testcross performance for GW. Moreover, the 
GW trait appeared to exhibit lower heterosis in testcross Fj 
population in both environments. The negative correlation 
between IL performance values and H MP values of the Fj 
testcross population clearly indicated that additive and dom- 
inant gene action acted independently in the testcross popu- 
lation. Furthermore, QTL and HL analyses demonstrated 
that six of the 41 HL were detected in QTL analysis and ex- 
hibited less genetic overlap with QTLs, consistent with re- 
sults reported by Hua et al. (2003); therefore, heterosis and 
trait performance may be conditioned by different sets of loci. 

A number of trait-enhancing QTL alleles have been re- 
ported; however, few favorable HL have been identified. Uti- 
lizing heterosis from the cross between indica and japonica 
is a trend in hybrid rice breeding programs. Due to the com- 
plexity of heterosis, conditioned by various factors includ- 
ing genetic background and environments, the detection of 
favorable HL between indica and japonica is a new chal- 
lenge for rice breeders. Our research indicates that introgres- 
sion lines are a good tool to identify HL between indica and 
japonica. ILs give a more reliable implication for genetic 
improvement, mainly by benefitting from more precise esti- 
mation of genetic effects against a relatively uniform and 
probably elite background (Tanksley and Nelson 1996). For 
the population used in this study, derived from a cross be- 
tween indica and japonica, hybrid weakness, sterility, and a 
large variation in heading date were observed. It might be 
possible to detect such recessive loci as heterotic loci. In 
fact, three loci {qGp5\ qSsp2b and qGwlb) appeared to 
have negative effects in IL analysis, but positive effects in 
H MP analysis. Furthermore, the qSsp2b locus was located in 
a position similar to the hybrid breakdown locus hbd2 re- 
ported by Matsubara et al. (2007). The three HL above are 
not viable candidates for the improvement of rice yield po- 
tential; however, analysis of the heterotic genetic basis of 
these recessive loci was viable. In this study, except for the 
three HL above, 24 HL showed significant positive heterotic 
effects on yield-related traits. In particular, grain yield per 
plant, the number of filled grains per panicle, and the num- 
ber of panicles per plant showed positive heterosis. Further 
study including fine mapping of these HL and cloning for 
yield-related traits will facilitate the use of these characteris- 
tics in breeding programs. 
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